function xx=fun1(t,x) 
a=1;b=0.1;c=0.5;d=0.02;
xx=[x(1)*(a-b*x(2));x(2)*(-c+d*x(1))];
